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ОСОБЕННОСТИ ВЫБОРА ПРЕДСТАВЛЕНИЙ ПОЛЕЙ й 
ВСПОМОГАТЕЛЬНЫХ ИСТОЧНИКОВ ЧАСТИЧНЫХ ОБЛАСТЕЙ 


В статье показана возможность уменьшения вычислительных затрат при решении граничных задач за счет 
рационального представления функций Г рина. На примере решения задачи электродинамического анализа 
проиллюстрирована ВОЗМОЖНОСТЬ оптимальной алг ебраизации математических соотношений, 

Ключевые слова: рупорный излучатель, продольно-слоистое заполнение, метод Г. алеркина, интегрирование. 


Введение. При решении граничных задач могут возникнуть вычислительные сложности, связан- 
ные с особенностями математического представления соответствующих функций, описывающих 
волновые процессы, протекающие в системе. 

Постановка задачи. На примере решения задачи электродинамического анализа ступенчатого 
рупора с продольно-слоистым диэлектрическим заполнением проиллюстрирована возможность 
уменьшения вычислительных затрат за счет рационального представления функции Грина (поля 
вспомогательного источника). 

Уменьшение вычислительных затрат при решении граничной задачи. В рассматриваемой 
модели излучателя (рисунок) слои диэлектрика (области \/1, \3, \/4) являются замкнутыми облас- 
тями, а свободное полупространство (область \2) и полубесконечный волновод (область \0) — 
полуоткрытыми областями. 

















Модель рупорного излучателя 


В соответствии с методом частичных областей [1] полагаем, что каждая частичная об- 
ласть электромагнитно изолирована, сторонние источники в областях отсутствуют, а возбуждение 
областей осуществляется эквивалентными токами, распределенными на поверхностях, ограничи- 
вающих области. Накладываем на тангенциальные компоненты электрического поля вспомога- 
тельного источника условие равенства нулю на поверхности, ограничивающей область. Тогда по- 
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ле в частичной области возбуждается эквивалентным магнитным током. Поскольку на идеально 
проводящих участках поверхности, ограничивающей область, касательная составляющая вектора 


Е равна нулю, эквивалентные магнитные токи, пропорциональные Е ‚ отличны от нуля только 
на неметаллизированных участках поверхности: 


БН (р) = [7"(@)Н“(р,9)а$,. (1) 


Здесь Н“” (р, 4) - возбуждаемые в точке р векторы напряженности магнитного поля вспомога- 
тельного магнитного источника (пропорциональные функции Грина данной области), распо- 


ложенного в точке а, находящейся на поверхности, ограничивающей область; И (а) — эк- 
вивалентный магнитный поверхностный ток. 


=, Е (а) (2) 
где и - внешняя нормаль к поверхности 5’, ограничивающей область; Е(а)- вектор напряжен- 


ности электрического поля в точке а; Ь - единичный вектор, вдоль которого ориентирован 
вспомогательный магнитный диполь. 

Полагаем, что излучатель возбуждается волнами типа Т, набегающими по волноводу (об- 
ласть \0) из минусбесконечности. В связи с этим в системе существуют только следующие со- 
ставляющие электромагнитного поля: Н,, Еу, Е.. 

Для примера запишем представление для составляющей Н, магнитного поля в области \1. 
Принимая во внимание совпадение смежных поверхностей 53, и 55, 54, И 58, 52, И би, а также 
совпадение неметаллизированного участка поверхности $! с 5, имеем в точке наблюдения 
(Ру) частичной области \1: 


Н,„ (ри) = | Е», (а»)Н и, (руль 4х 95% — 
50 


-. | Е, (Ч5, Н (Ри ,Ч5, )45> их. 


52 


№ | (45, НН. (Риь95, 7455 + 
56 
+ [Е 


2у4 
58 


(45, )Н (Руль 95, 4$. (3) 


Здесь Ё Е Е РЕ - тангенциальные составляющие электрического поля об- 
Уто Ур! уз уд 


ластей \0,№1,\3 и \4 в точках и поверхностей 55, 52, 55 и 58, ограничи- 
Е 


вающих область \1. 

Учитывая, что каждую поверхность можно рассматривать как плоскость, перпендикуляр- 
ную продольной оси волновода (направленного по осям у либо 2), составляющие электрического 
поля целесообразно представлять в виде поперечных компонент электрического поля соответст- 
венно ориентированного волновода [1, 2]. 

Аналогичным образом можно представить поля в областях \2, \З и \4. 

Поле в области \0 представлялось в виде разложения по системе волноводных гармоник. 

Для обеспечения единственности решения электродинамической задачи необходимо при- 
равнять тангенциальные составляющие магнитного поля смежных областей на границах их раз- 
дела, т.е. совместить точки наблюдения на границах раздела смежных областей. В результате 
получим систему интегральных уравнений, которая в операторном виде представлена следующим 
образом: 
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Н„, 02 =0)=Н,, (7,2 =0), 

Н„,„ (2 =0)=Н,. 0,2 =0), 

Н,, (0,7 =0)=Н,,(у,2' =0), 

НО =6:,2') =Н,„, (у = 6:,2'), (4) 
НО’ =В+,2)=Н,„.(у=Ь+Ь,2)), 

Н „0, =-4)=Н„„ (у, =-4), 


где штрихованные у’ и 2’ обозначают координаты точек наблюдения. 





Таким образом, приравниваются как тангенциальные магнитные, так и тангенциальные 
электрические компоненты полей смежных областей. 
Решение системы интегральных уравнений выполнялось методом Галеркина. 


К м 
Для конкретизации представления полей Н, (р,4 ) в частичных областях необходимо 


решить задачу о возбуждении соответствующей частичной области бесконечной нитью магнитно- 
го стороннего тока, ориентированной вдоль оси Ох и имеющей единичную амплитуду и нулевую 
фазу. Решение производится в два этапа. 

На первом этапе решается задача о возбуждении плоского волновода. Решение неодно- 
родных уравнений Максвелла сводится к решению неоднородного уравнения Гельмгольца для 


комплексных амплитуд магнитного векторного потенциала А” [2]. Вектор напряженности маг- 
нитного поля выражается через векторный потенциал. Решение уравнения Гельмгольца пред- 
ставлялось в виде разложения в ряд Фурье для волны, бегущей в направлении оси 2 (первый 
волновод) и оси у (второй волновод): 


А» (9,2) = Увы 00$ (рлЬ у) с0з (рлЬ"у’ехр [+у.„(2-2) |, (5) 
1 р=0 

4, (,2) = Увьупьсо (ря 4 "=)с0з (рла“=')ехр[+у»(’-5)]. = (6) 
р=0 


Здесь: 






пр р=0 
пр р>=0. 
Нештрихованные у и 2 соответствуют координатам источников. 
На втором этапе к полученному выше решению неоднородного уравнения добавляется 


решение однородного уравнения Гельмгольца. Это решение для каждой парциальной волны со- 
держит два неизвестных коэффициента : 


м 1 В = ЙИр2 1р2 — ей 
А’одн 25. У = (Еве 77+ Е, е”””) сз. ( рл В, 1у)соз-(рлЬ у’); (7) 
1 р=0 
1х а о 
Агодн = аа (Еы,е и" + Е. „е“””)со$: (рл 4 '2)с03: (рл 417’). (8) 
р=0 


Неизвестные коэффициенты Е определяются из граничных условий для тангенциальных 
составляющих напряженности электрического поля на торцевых стенках резонатора. В результа- 
те получим следующие выражения для магнитного поля вспомогательного источника резонатора 
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(область\1), образованного из волновода, в котором распространяется волна в направлении оси 2 
или у (резонатор с продольным направлением по оси гили у): 


‚1х 8 рл рл 
Н" =—= Г с0$ (у-Ь о > 
й . В >. 7.5й(У А) | В, | В, } 





сй(у, „2)сй [5 „(а +2) | при 2' <2 











. (9) 
сй(у,,7')сй [у.„(а + 2) | при 2.> 2; 
е — 4) 
Н*" = 1 у =, (-1) соь| 27 2 2 =], 
т Ч ро 1 ‚51 ,,В,) а а 
сй(уг,(у-В))сй а (В в) | приу'>у (10) 


сй(у, ‚(Вей у, (-Ь-Ь,) | приу" < у. 


Таким образом, соотношения для полей вспомогательных источников могут быть пред- 
ставлены различными выражениями. 

Выбор выражения для поля вспомогательного источника целесообразно увязывать с ви- 
дом представления эквивалентного тока, распределённого на конкретной поверхности, охваты- 
вающей область. Рассмотрев два варианта представления полей вспомогательных источников, 
установили, что более предпочтительным является такое представление поля вспомогательного 
источника, которое характерно для резонатора с продольным направлением, совпадающим с 
продольным направлением волновода, через поперечные составляющие электрического поля ко- 
торого выражается эквивалентный поверхностный ток, входящий в подынтегральные выражения. 

Проведем вывод коэффициентов, описывающих связь гармоник поверхностного тока на 
5, с гармониками поверхностного тока, распределенными на поверхности $3, при сшивании по- 
лей областей \/0 и \1 по границе раздела $, Эта группа коэффициентов вычисляется после соот- 
ветствующих преобразований третьего слагаемого правой части выражения (3). При применении 
метода Галеркина в качестве весовых функций использованы базисные функции разложения по- 
тул 
Ь 





верхностных токов на поверхности 5: с ’->) ‚ базисные функции на поверхности $3 
имеют вид с05| 7 - 2| где Тзг - номера гармоник волновода, ориентированного вдоль 


, . 
оси у, а И - номер весовой гармоники. В соответствии с методом Галеркина после умножения 


интегрального уравнения, определяющего связь поля возбуждающего волновода (область \0) с 
полем в области \1 (см. шестое уравнение системы (4)) на весовые функции, производится ин- 
тегрирование по точкам наблюдения. 

В первом случае при выводе коэффициентов используем выражение для поля вспомога- 
тельного источника резонатора с продольным направлением по оси 2. Точки наблюдения при 


этом имеют координаты у, 7= - 4 точки источников — У = (3,2). Подставляем выражения для 


электрического поля и поля вспомогательного источника при условии, что 2’<2. В результате по- 
лучим следующее выражение: 
‚А ту л ее т. 
т 0 ' ЗГ 
С10Г, 65 | С05 В (у = у.) | (-—)>, УЗгиь с05 а г} 
о 4 


У 3 тзг 





- (11) 


1 = рл 
х(—1=,) - с0$ (-Ь,) | сй(у, „2)а2ау’. 
СВ рр И я 
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После вычисления интегралов получим: 








Я 5 
ОЕ Е 
г =3Б, р= бу? + тзгл (12) 
ы а 
где 
рл 
р т. рт : рл 
| и :| 1) |2 (ву, ь) “х (5. м) 
ЕР . у 
| В и 
в В 
Е1 = 
—2 с0$ р, — у) при РА = 20, 
, в В 
рл _тл 
В при ==“ = 
. В Ь 





Во втором случае используем выражение для поля вспомогательного источника резонато- 
тзг В 


т 
ра с продольным направлением по оси у. Исходное выражение для коэффициента С101 


= ») Г лы |" и. 


163 тзг =0 


этом случае запишется: 


, и +5 
С101^ = | к" 











я —4 
И ЕН [2 (13) 
- й 5. -ЁЬ, } | 42а". 
^( Е а ре [Ин-( 3 ‚) — 
После интегрирования: 
| 1 
С (1х 
31 а ЗГтзг 5й( БАЙ | "=" 
17тз г 
о (14) 

о [Узлы (4+5 -В, В) |--5® [уши (и -В-5, УЕ, 


1 при и -0: 
1 
в - тзг = 0: 


1 =0 
© = 
” [2 при т. #0; 


где 
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Аналитические выражения для вычисления коэффициента ©1091. различны. Соотно- 
шение (12) включает в себя суммирование в бесконечном пределе по индексу р, в то время как 
соотношение (14) позволяет вычислять значение коэффициентов непосредственно. Однако чис- 
ленные значения коэффициентов, вычисленные по соотношению (12), при увеличении числа 
членов суммирования по индексу р стремятся к численному значению коэффициентов, вычислен- 
ных по (14). 


, Й 
В таблице приведены значения коэффициентов С10Г, при различных то и 3г, вы- 


численные по соотношению (12) при учёте разного количества слагаемых ряда по р. В последней 
справа колонке таблицы приведены значения коэффициентов, рассчитанные по соотношению 
(14). Расчеты проводились при следующих параметрах излучающей системы (см.рисунок): 


В =0,30л, В, =В:=6=0,23п, 4=0.8п, &=1, а=3, &=1, в=а=3, У] = 0.195л. 


Значения коэффициентов ©1091”, 

















Р 10 50 100 150 500 С10 [то 

тзг 

Сто =0 0,00000000- 0,00000000- 0,00000000- 0,00000000- 0,00000000- 0,00000000- 
тзг=0 -0,538781237]| -0,539275190] -0,539281207] | -0,539281753] | -0,539281752] -0,539281752] 
т'=0 0,233486849+| 0,232470465+ 0,232457916+ | 0,232456777+ | 0,232456780+ 0,232456780+ 
тзг=2 +0,00000000] | +0,00000000} +0,00000000] +0,00000000] +0,00000000] +0,00000000] 
т =1 0,00000000+ 0,00000000+ 0,00000000+ 0,00000000+ 0,00000000+ 0,00000000+ 
тзг =0 +0,236506601| +0,236143398} | +0,236137785] | +0,236137316] | +0,236137353] | +0,236137353} 

] 

т =1 0,129104918+| 0,1283556224-+ | 0,128344520+ | 0,128343542- | 0,128343619+ 0,128343619+ 
тз-=2 +0,00000000] | +0,00000000} +0,00000000] +0,00000000] +0,00000000] +0,00000000] 





























Как следует из таблицы, при увеличении числа членов ряда по р в выражении (12) значе- 


т'о 
ния коэффициента С10Е,‚. стремятся к значениям, вычисленным по выражению (14). Для вы- 


бранных размеров излучающей системы совпадение значений коэффициентов в пятом знаке по- 
сле запятой имеет место при р>100. При увеличении размера Ё: совпадение значений с заданной 
точностью наступает при учёте меньшего количества слагаемых в выражении (12). 


Выводы. Результаты вычисления подобных элементов блочной матрицы показали, что числен- 
ные значения коэффициентов, рассчитанные по соотношениям (12) и (14), совпадают. При этом 
совершенно очевидно снижение вычислительных затрат в случае использования поля вспомога- 
тельного источника резонатора с продольным направлением, совпадающим с продольным на- 
правлением волновода, по поперечному размеру которого производится интегрирование. Это 
обусловлено тем, что хотя соотношения для расчета коэффициентов имеют один порядок слож- 
ности, но в одном случае (см. формулу (14)) вычисления производятся один раз, а в другом слу- 
чае (см. формулу (12)) — сотни раз. 
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